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I. ii^TPfinurTlON 


The Goddaix) Laboratory for Atmospheres (GLA) has over the past few > 
acquired a large data base of OOZ and 12Z initialized analyses (1980-87) from the 
Center for Medium Range Weather Forecasts (ECMWF). In support of our on-; 
research activities we have found it convenient to compile these analyses in the 


fltias 

The original ECMWF/WMO data was obtained on a 2.5 ° latitude by 2.5 
grid in a packed format. This data was unpacked and interpolated to a 4 0 latitud. 
longitude grid for the purpose of computing various statistics and plotting. App 
describes furthur details of the data and data processing. In part I of this atlas 
various first moment quantities for monthly, seasonal and annual averages on a 
cylindrical projection as well as cross section maps of zonal averages. We also 
average seasonal cycle and anomalies during the 8 year period for selected quar 
II presents various second moment quantities including both transient and stand 
components. This includes some time filtered quantities which separate the con 
of the band pass (2.5 - 6 days) and low pass (> 10 days) regions of the spectre 
The ECMWF analysis system has undergone a number of changes dun 
period 1980-87. Trenberth and Olson (1988) have done an extensive study of . 
analyses for 1979-86 and examined the impact of analysis system changes on o 
(see their Table I and discussion). The change to diabatic nonlinear normal me 
initialization in September of 1982 dramatically increased the strength of the tro 
divergence. For quantities sensitive to the divergence we have limited the time 
1984-87 For these fields the shorter time-averaging period is clearly indicates 
the averages include all 8 years. It must be emphasized that many other change 
to the analysis system during this time period leading to inhomogeneities in the 
is particularly true for the moisture fields (see eg. Trennberth and Olson, 198,' 
attempt has been made to take this into account when performing the time aver 

We note that other recent compilations of assimilated data that may be 
comparisons are those of Trenberth and Olson (1988) based on essentially the 
(1979-86) and Lau (1984) which examines the FGGE year based on analyses 

GFDL. 

In accordance with ECMWF conditions any publications making use ( 
and/or data in this atlas must acknowledge ECMWF as the data source. 


V 


Btt fc M / 


;ars 

European 

aing 

rm of an 

longitude 
by 5° 
ndix A 
: present 
lobal 
resent the 
ites. Part 
ng eddy 
ributions 
a. 

g the 

ie ECMWF 
ita quality, 
e 

>ical 

averages to 
Otherwise 
> occurred 
data. This 
). No 
ges. 

seful for 
ame data 
iroduced at 

* the figures 


JNIENTlQNAbU 1 i 


PRECEDING PAGE BLANK NOT FILMED 


II. LIST OF pimYTfTipfi 


The follow, ng l.s,s the quantities and units used for plotting. The contour intervals 

”“' y , frnm thC P ' 0,S ,hemselves - The global maps consist primarily of the 850 
and 200 mb Helds. Exceptions to this are the heigh. Held which is presented a. 300 mb the 

7™ 7'°7 y Whlch 15 presen,ed 31 500 “> d 'he moisture field which is shown only 

* T e cros s s e cti°n s are based on all the levels obtained from ECMWF which are 

’ 50, 700 ’ 50 °’ 300> 200, and 100 mb - P,ots ^ sh own for monthly, seasonal and 
annual averages. Our organization of the fields follows essentially that of Lau (1984) 

The basic quantities and units (for the plots) are 


1) U 

2) V 

3) T 

4) Z 

5) Ze 

6) Q 

7) to 

8) y 

9) X 

10 ) 


east/west component of wind 
north/south component of wind 
temperature 
geopotential height 

eddy geopotential height (zonal mean removed) 

specific humidity 1 

vertical velocity 

stream function 

velocity potential 

sea level pressure 2 


m/s 

m/s 

°K 

m 

m 

g/kg 

nanobars/s 
m 2 /s x 10 6 
m 2 /s x 10 6 
mb 


With- Th t flddS “7 gr0UPed aCCO,din8 10 Whe,her ,hey 3,1 m ° n,hly - or annual averages 

bm tetsT H f r maPS PreSe " ,ed f,rS ‘ “ ^ 2 °" al "»•" * 

([. .]). The final maps present the seasonal cycle and anomalies (monthly deviations from 

* cycle) of sclected zo " al avera 8 e d fields. We indicate below the quantities and (for 
convemence) the relative positions of the plots in the atlas for each data group 


'The original field was relative humidity. See Appendix A for the 
See Appendix A for the sea level 


conversion to specific humidity. 


pressure calculation, from the height and temperature field. 




The following quantities describe the seasonal cycle which is the average O' tr all years 
(1980-87). We also show the difference between the first and last 4 years . .984/87 
1980/83). These quantities appear after the annual averages. 



vu 




The anomalies are the monthly averaged deviations from the seasonal cycle. For [V] 
[w], and [Q] we remove separately the seasonal cycle for 1980-83 and 1984-87. 


[U]* 

Quantify 

200 mb anomalies 

Relative 

1 

2 

[U]* 

850 mb anomalies 

[V]* 

200 mb anomalies 


[V]* 

850 mb anomalies 

4 

[CO]* 

500 mb anomalies 

5 

[TJ* 

200 mb anomalies 

mJ 

6 

[T]* 

850 mb anomalies 

7 

IQ}* 

850 mb anomalies 

8 


Appendix B describes how to access these data sets on the Amdahl V7 and V6 
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IV. APPENDIX A- THE DATA PROCESSING 


The ECMWF analyses are the initialized fields provided to GLA (code 61 1) in an 
internal bit format. The data processing involves several steps including unpacking and 
converting to a 2 X 2.5 degree grid. During this step the relative humidity is converted to 
specific humidity and the sea level pressure is computed. Then an interpolatiion is made to 
a 4 X 5 degree grid. The stream function and velocity potential are computed from the 
lower resolution monthly-averaged fields. Both resolutions of the basic parameters are 
available on tape but all plotting is done from the lower resolution fields. Further details of 
the data are provided below. 

ORIGINAL PACKED DATA AS RECEIVED FROM ECMWF: 

2.5° latitude by 2.5° longitude in internal bit format 
1000, 850, 700, 500, 300, 200, 100 hPa 
separate 00Z and 12Z tapes 


DATA (IM,JNP) 

IM= 1, .. , 144 starting at 177.5 W and going east to 180 E 
JNP = 1, .. , 73 starting at North Pole 

(1) geopotential height (m) 

(2) temperature (°K) 

(3) u (m/s) 

(4) v (m/s) 

(5) relative humidity (%) 

(6) vertical velocity (Pa/s) 


UNPACKED ECMWF DATA (Version 8 format) 

2° latitude by 2.5° longitude or 4° latitude by 5 0 longitude 
In version 8 format (see documentation- VERSION 8 I/O utilities) 
1000, 850, 700, 500, 300, 200, 100 hPa 
00Z and 12Z combined for each tape 


DATA(IMJNP) 

IM=1, .. , 144 (72) starting at 180° W and going east to 177.5° E (175° E) 
JNP= 1, .. , 91 (46) starting at South Pole 


x 



Quantity 


(1) sea level pressure (hPa) 

(2) geopotential height (m) 

(3) temperature (°K) 

(4) u (m/s) 

(5) v (m/s) 

(6) specific humidity (g/g) 

(7) vertical velocity (Pa/s) 

(8) divergence 


The specific humidity (q) is computed from the relative humidity 
q = q s * RH * 0.01 


where 

q s = 0.622 * e s / (p - .378 e s ) 

and e s is the saturation vapor pressure determined from the empirical formulae 

6.1078 * EXP { 21.8747 * ( T- 273. 16)/ (T - 7.66)} forT<273. 6 

6.1078 * EXP { 17.269 *( T- 273.16 )/ (T - 35.86) } forT> 273 16 


The computation of sea level pressure (p s i) is based on the hypsometric equation nd an 
assumed lapse rate (see eg. Trenberth , JGR, 86, C6) so that 

p sl = 1000* EXP (9.81* Z/R T v ) 


where Z is the height of the 1000 hPa surface and the mean virtual tempe iture (T v ) 
is given by 

T v = T * (1 + 0.61 * q ) + 0.5 * (0.0065 * Z ) 


xi 



V. APPENDIX B- SOURCES OF DIGITIZED DATA 


The data fields (for Part I of this Atlas) are currently stored in the MVS partitioned 
data set (PDS), W3CKP ECMWF ATLAS and are meant for use by NASA scientists with 
access to the Amdahl mainframe. This file will also be made available on the IBM 3081 . 
We also anticipate that these fields will be made available on floppy disks for both IBM 
compatible and Macintosh personal computers in the near future. For further information 
concerning the use of this data contact 

Siegfried Schubert 
NASA/GSFC Code 61 1 
Greenbelt, MD 2077 1 
(301) 286-3441 


1 . MVS JCL and CMS EXEC 

a) //FT11F001 DD DSN =W3CKP.ECM WF. ATLAS(xxxxxx),DISP=SHR 

b) FI 11 * DSN W3CKP ECMWF ATLAS (MEMBER xxxxxx 

Here, xxxxx is the PDS member name. 

2 . Read format and PDS member names 

Longitude/Latitude Fields 
MEMBERS: 

SLP sea level pressure (mb) 

U200 u-wind at 200 mb (m/s) 

U850 u-wind at 850 mb (m/s) 

V200 v-wind at 200 mb (m/s) 

V850 v-wind at 850 mb (m/s) 

Z300 geopotential height at 300 mb (m) 

T200 temperature at 200 mb (°K) 

T850 temperature at 850 mb (°K) 

Q850 specific humidity at 850 mb (g/g) 

W500 omega at 500 mb (Pa/s) 

VEL200 velocity potential at 200 mb (m 2 /s) 

VEL850 velocity potential at 200 mb (m 2 /s) 

STF850 stream function at 850 mb (m 2 /s) 

STF200 stream function at 200 mb (m 2 /s) 

REAL FLD(72,46) 

DO 100 IMON = 1, 12 
100 READ(n) I, FLD : MONTHLY MEAN (JAN-DEC) 

DO 200 ISN =1,4 

200 READ(n) INDEX, FLD : SEASONAL MEAN (DJF-SON) 

READ(n) INDEX, FLD : ANNUAL MEAN 

FLD(I,J) J= 1 is at the South Pole and 1=1 is at 180° W (increasing eastward to 175 E) 


(1984-87 average) 
(1984-87 average) 
(1984-87 average) 


xu 



Latitude/Pressure Sections 
MEMBERS: 


LTPQ 

LTPT 

LTPU 

LTPV 

LTPW1 

LTPW2 


zonally-averaged specific humidity (g/g) (1980-1 7) 
zonally-averaged temperature (°K) (1980-87) 
zonally-averaged u-wind (m/s) (1980-87) 
zonally-averaged v-wind (m/s) (1984-87) 
zonally-averaged vertical velocity (Pa/s) (1980-83 
zonally-averaged vertical velocity (Pa/s) (1984-87 


DIMENSION LTP(46,7) 
DO 100 IMON = 1, 12 
100 READ(n) INDEX, LTP 
DO 200 ISN =1,4 
200 READ(n) INDEX, LTP 
READ(n) INDEX, LTP 


: MONTHLY MEA I (JAN-DEC) 

: SEASONAL ME/ V (DJF-SON) 
: ANNUAL MEAN 


LTP(J,L) J= 1 is at the South Pole and L=1 is at 1000 mb. The other pr< .sure levels are 
850mb, 700mb, 500mb, 300mb, 200mb, lOOmb. 

Seasonal cycle climatolog y 

MEMBERS : CLI8087, CLI8083, CLI8487 


DIMENSION CLI( 12,46) 

DO 100 IV = 1, 8 
READ(n) INDEX, CLI 
100 CONTINUE 

CLI(K,J) J= 1 is at the South Pole and K=1 is January. 

Quantity 

zonally-averaged u-wind at 200 mb (m/s) 
zonally-averaged u-wind at 850 mb (m/s) 
zonally-averaged v-wind at 200 mb (m/s) 
zonally-averaged v-wind at 850 mb (m/s) 
zonally-averaged temperature at 200 mb (°1 ) 
zonally-averaged temperature at 850 mb (°1 ) 
zonally-averaged specific humidity (g/g) 
zonally-averaged vertical velocity (Pa/s) 


1Y 

1 

2 

3 

4 

5 

6 

7 

8 


i 


xm 


mi’IHIh 


Departure from seasonal cycle 

MEMBER : DEP8087 

DIMENSION DEP(96,46) 

DO 100 IV = 1, 8 
READ(n) INDEX, DEP 
100 CONTINUE 

DEP(K,J) J= 1 is at the South Pole and K=1 is January, 1980. 


IX Quantity 


1 

2 

3 

4 

5 

6 

7 

8 


zonally-averaged u-wind at 200 mb (m/s) 
zonally-averaged u-wind at 850 mb (m/s) 
zonally-averaged v-wind* at 200 mb (m/s) 
zonally-averaged v-wind* at 850 mb (m/s) 
zonally-averaged temperature at 200 mb (°K) 
zonally-averaged temperature at 850 mb (°K) 
zonally-averaged specific humidity* (g/g) 

zonally-averaged vertical velocity* (Pa/s) 


* These anomaly fields are based on separate climatologies for 1980-83 and 1984-87. 
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VI. THE FIRST MOMENT FIELDS 
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fEC^r^M ™ atl f ° f f urope ^ n Center for Medium Range Weather Forecasts 
(ECMWF) initialized analyses for 1980-87. Various first moment quantities are 

presented for monthly, seasonal, and annual averages on a global cylindrical 
projection, as well as cross section maps of zonal averages. Global maps of 

S ^ eam functlon ’ a nd velocity potential are presented at 
850 and 200 mb. In addition, we include global maps of the 300-mb height field 
(total and eddy), the 500-mb vertical velocity, the 850-mb moistura field and sei 
level pressure. The average seasonal cycle and anomalies during the 8-year peri ^ 
are presented for selected quantities. * V 
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